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Abstract. Perpetual American options are financial instruments that can be readily 
exercised and do not mature. In this paper we study in detail the problem of pricing 
this kind of derivatives, for the most popular flavour, within a framework in which 
some of the properties — volatility and dividend policy — of the underlying stock can 
change at a random instant of time, but in such a way that we can forecast their 
final values. Under this assumption we can model actual market conditions because 
most relevant facts usually entail sharp predictable consequences. The effect of this 
potential risk on perpetual American vanilla options is remarkable: the very equation 
that will determine the fair price depends on the solution to be found. Sound results 
are found under the optics both of finance and physics. In particular, a parallelism 
among the overall outcome of this problem and a phase transition is established. 



PACS numbers: 89.65. Gh, 02.50.Ey, 05.40.Jc, 02.30.Jr 
1. Introduction 

Pricing financial derivatives is a main subject in mathematical finance with multiple 
implications in physics. In 1900, five years before Einstein's classic paper, Bachelier [1] 
proposed the arithmetic Brownian motion for the dynamical evolution of stock prices 
with the aim of obtaining a formula for option valuation. Samuelson [2] noticed the 
failure of such market model that allowed negative values in stock prices, and introduced 
the geometric Brownian motion which corrects this unwanted feature. Within his log- 
normal model, Samuelson obtained the fair price for perpetual options, although he was 
unable to find a general solution for expiring contracts. The answer to this question 
must wait until the publication of the works of Black and Scholes [3], and Merton [I]. 
The celebrated Black-Scholes-Merton formula has been broadly used by practitioners 
since then, mainly due to its unambiguous interpretation and mathematical simplicity. 

This mathematical simplicity in the Black-Scholes-Merton scheme has nevertheless 
a drawback: the model poorly adapts to those evolving conditions that affect actual 
derivatives and real markets present. In particular, empirical analyses conclude that 
volatility, roughly the diffusion parameter, must be considered as a changing (random) 
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magnitude rather than a mere constant, as in the Black-Scholes-Merton formula. Many 
models have been developed in this direction, but among them a few deserve special 
emphasis because their historical imprint: the works of Hull and White [5], Wiggins [6], 
Scott [7J, Stein and Stein [8|, or Heston [9] belong to this selected group. They 
are collectively termed as stochastic volatility models because volatility becomes a 
continuous process that follows its own stochastic differential equation. Quite in the 
opposite direction we find another possible approach to the issue: the Markov-modulated 
geometrical Brownian motion models. Within these models the market coefficients 
change in a deterministic way but at random times, and according to that they are 
generically known as regime- switching models. Behind these lies the theory that the 
value of the parameters depends on the state of the economy, that suffers from seasonal 
changes. Seeds of this idea appeared in Barone-Adesi and Whaley [10], but the first 
time a model of this kind was properly settled is in the work of Naik [llj. Since 
then such models have been used to discuss European [121 H31 EEU US], Russian [16] or 
American [T7[ [TBI [13 EO, [21] option prices, but also in asset allocation and portfolio 
optimisation problems [22[ [25[ l2T[ 1251 126] . 

Once we have mentioned different option flavours, we must also recall that the 
Black-Scholes-Merton formula is only is applicable to European options — derivatives 
that can be exercised at maturity alone — , whereas most of the exchange-traded options 
are American — they can be exercised anytime during the life of the contract. Kim [27] 
provided an integral representation of the American vanilla option price, but the explicit 
solution remains still unknown. 

The real difficulty behind pricing American options lies in the fact that one has 
to solve free boundary problems for partial differential equations, equations which 
sometimes have clear physical interpretation: McKean [28], in what is the earliest 
approach to the issue, unambiguously speaks about the heat equation in an economic 
publication. Here the boundary represents the optimal exercise price, the stock price 
that triggers the early exercise of the option, and in general one can determine it just once 
the pricing expression is known, what leads to a circular problem. As a consequence, only 
restricted circumstances allow for closed formulas in American-like problems, whereas 
in the most general scenario analytical or numerical approximate methods must be used 
instead. One of theses favourable instances corresponds to a major simplification in 
the American problem: the assumption that the option never expires. Non-expiring 
or perpetual options serve as a good stating point in the resolution of the complete 
problem [10] because the absence of maturity usually removes any explicit temporal 
dependence that those differential equations present. 

In this article we are going to tackle the problem of pricing perpetual American 
vanilla calls and puts within a market model that can be considered as a degenerated 
instance of a regime-switching model with only two states and where one of the transition 
probabilities is set to zero: We will let volatility and dividend rate perform a single 
regime change at some unknown instant in the future. The value before and after the 
regime change of such magnitudes are assumed to be known in advance, what indeed 
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determines if the change in the dynamics has taken place or not. 

Part of the present analysis is formally included as a limiting case in the paper by 
Guo and Zhang cited above [H] , where the authors studied perpetual American vanilla 
put options with regime switching in absence of dividends, but there are divergences 
as well. On the one hand, we will also pay attention to the call case. Since dividend- 
paying assets are considered, a perpetual call may differ in price from its underlying, even 
though the option must not be early exercised. On the other hand, their study on puts is 
complete and very general but few explicit details about the closed-form solution found 
are given. This is specially relevant in the analysis of the properties of the equations that 
determine the optimal exercise price: In essence one obtains different pricing formulas, 
each of them coming from an extreme value problem, and must elucidate the right choice 
at every moment. In spite of the multiplicity of parameters that our model involves, we 
will show how the value of a single magnitude answers the dilemma. Therefore, there is 
a formal resemblance between this behaviour and a phase transition in thermodynamics. 

The paper is structured as follows: In Section [2] we present the market model, 
the securities traded in it and their general properties. In Section [3] we introduce 
the concept of hedging portfolio and show how it can be used for pricing derivatives. 
Section @] is specifically devoted to perpetual American vanilla options, we settle the 
properties of these ideal derivatives, and stress their interest from the point of view of 
finance and physics as well. In Section [5] we quote the analytic expressions found, and 
emphasize the interpretation of the optimal exercise price as the outcome of an extremal 
problem. In Section [6] we discuss the most appealing properties and particularities of 
the formulas that conform the solution to the problem, by illustrating our inferences 
through graphical examples. Conclusions are drawn in Section 0, and the paper ends 
with a lengthy Appendix where we present step-by-step calculations for call options. 

2. The market 

Let us begin with the general description of the stochastic properties of the securities 
that are traded in the market. The first security to be considered is a zero-coupon bond, 
a riskless monetary asset with a market price B t with deterministic evolution: 



where r, the risk-free interest rate, is assumed to be constant and positive. 

The second security present in the market is the stock S t , whose time evolution for 
t > fulfils the following Ito stochastic differential equation (SDE): 



where W t is a Wiener process, a one dimensional Brownian motion with zero mean and 
variance equal to t. The drift, \i ty and the volatility, Ot > 0, are stochastic processes as 
well: we will assume that their initial values are fi a and o~ a , and after that moment they 



B t = B e r \ 



(1) 



dS t = fit-St-dt + at-St-dWt 



(2) 
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may simultaneously change to some other (different) fixed values, /if, and o^, known in 
advance. Such a turnover can take place only once in a lifetime: 

/i t = /i a l i<T + /ifel^r = Ha + - A«a) !t>r, (3) 

0~t = 0- a l t<T + Ofclt>r = CTa + (0"6 - CTa) ltJST- (4) 

We state that these magnitudes are stochastic because the instant r > in which 
the regime change occurs is a random variable. We will consider that r follows an 
exponential law with E p [r] = A -1 > 0. In mathematical terminology both processes are 
cadlag (right-continuous with left limits) semimartingales [29J E0] • In spite of the fact 
that fi t and at are not continuous processes, the price itself is continuous because we are 
not considering any jump contribution in (j2J). Therefore we may safely replace St- by 
St hereafter. Let us now analyse the meaning of (JHJ) and (J4j). The case of at is somewhat 
simpler since when a a ^ ab we may speak of a (elementary) stochastic volatility market 
model with just two possible volatility levels. The meaning of the change in fi t admits 
several interpretations [HI HHJ [25] , but we will concentrate our attention in the existence 
of two different (continuous-time) dividend pay-off regimes, fj, t — jit — $t, 

5 t = 5 a l t<T + 5 b l t ^ T = 5 a + (5 b - 5 a ) l t ^r, (5) 

with 5 a> b ^ 0. 

A third kind of securities are negotiated in the market: the derivatives of S. 
Derivatives [3T] are financial instruments whose value depends on (derives from) present 
and past quotes of one or more underlying assets, which commonly are simply referred 
as the underlyings. In particular, we are interested in options: contracts between two 
parties that give the buyer the right, but not the obligation, to buy (call) or sell (put) 
shares of the underlying stock at some prearranged price, the strike price K > 0, 
at or before some specific instant in the future, the maturity or expiration time T. 
More precisely, we will assume that the option depends on the underlying by means 
of the present asset value alone, S t . However, the option price does also depend on 
the stochastic process I t , It = lt> T , through the market parameters fit, a t and 5 t . The 
statement that an indicator function with a random variable in its argument constitutes 
a stochastic process may seem a delicate issue but it is mathematically correct: In fact, It 
is a submartingale under our measure and, by virtue of the Doob-Meyer decomposition 
theorem [29], it can be expressed as a sum of two terms, I t = A t + M t , the increasing 
adapted process A t , 

A t = A min(r, t), 

and the cadlag (local) martingale M t , 

M t = l- AE p [r|^]. 

In the previous expression Tt represents all the available information up to time t, 
and corresponds to what is known in mathematical terminology as a filtration: an 

J Throughout the text will denote the indicator function, which assigns the value 1 to a true 
statement, and the value to a false statement. 
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increasing one-parameter family of sub cr-algebras of T , Tt Q T s C T , t ^ s, where 
T is a cr-algebra of subsets of the sample space fl The sample space, T and the 
measure P — which appeared beside the expectation symbol — define the probability 
space (n,^ 7 , P). Finally, note that maturity introduces explicit temporal dependency in 
the value of the option as well, so we will have that P t = P(t, St, It)- 

In the present case in which the price of the option is a function of two stochastic 
processes, one of them with discontinuous sample paths, the appropriate Ito formula [30] 
reads: 

P(t,S t ,I t )=P a (0,S ) + [ t d s P(s,S s ,I s _)ds+ f d s P(s,S s ,I s _)dS s 

Jo Jo 

+ ^-^SSP^ ^ J -) ds + AP r (6) 

where 

AP T = P b {r,S T )-P a {r,S T ), 

and Pb(t, S t ) = P(t, S t , 1), P a (t, S t ) = P(t, S t , 0). We can infer the discontinuous nature 
of the option price from the presence of the last term in ([6]). The differential form of 
this expression is: 

dP t = d t Pt-dt + dsPt-dS t + l -a 2 t _S 2 t d 2 ss P t ^dt + AP t _dI t , (7) 

where the property P t - = Pit—, S t -, It-) = P(t, S t , J t _) has been systematically used. 
Depending on the notation used, one can find that the AP t _dI t term from the above 
expression is simply replaced with AP t - , but then one must define the random measure 
associated to the jumping process. We have decided not to proceed in that way because 
the appearance of dl t in the forthcoming expressions stresses the option price dependence 
on a stochastic process different from dS t . Therefore, it is a source of risk that cannot 
be explained in terms of the random evolution of the underlying asset. When part of 
the risk is not directly traded in the market, the market may be incomplete: we will not 
be able to reproduce the behaviour of the option by means of a replicating or hedging 
portfolio. The immediate consequence of such eventuality is that the risk premium 
coming from r is arbitrary to a certain extent, because investors can evaluate it on the 
basis of their own perceptions. We will see all of this in more detail in the next Section. 

3. Option pricing 

In this Section we will obtain pricing expressions for the derivatives of S under the 
market conditions specified above, by resorting to the hedging-portfolio technique. This 
approach is based on the idea that the fair price of an option must be equal to the 
value of some portfolio II made of different securities that mimics the behaviour of the 
derivative and hedges all the risk. If such a portfolio can be defined, the market model 
does not admit arbitrage [3l] . 

The first security to be included in the portfolio is the underlying asset, which will 
reproduce changes in the option price due to the evolution of the stock price S. The 
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second security in our portfolio is B: a long position in this security will provide a 
secure resort where to keep the benefits of an effective trading strategy, whereas a short 
position in these bonds will allow us to borrow money when we need it. These two 
securities cannot counterbalance all the stochastic behaviour of the price of the option 
as we have pointed out before: not all the influence of 5t and o t on the option price may 
be explained through St- Therefore, we need another security that can account for this 
contribution to the global risk. Since in the most usual situation markets do not trade 
such assets, we will consider the inclusion of a secondary option Q in the portfolio: a 
derivative of the same nature of P, but with a different set of contract specifications. 
In particular we will assume that Q has a striking price K', different from K. This is a 
standard technique [32j [33] intended to complete the market. 

Summing up, our aim is to reproduce the evolution of Pt by means of a portfolio 
n t built out of a mixture of u t shares S t , <fit units of the riskless security B t , and ipt 
secondary options with value Q t : 

lit = v t S t + <p t B t + AQt, (8) 

in such a way that IT = P t holds almost surely. The predictable caglad (left-continuous 
with right limits) processes u t , <p t and ip t constitute what is known in mathematical 
finance terminology as a trading strategy, and the capital gain T> t associated to it is 
given by: 

V t = [ u t dS t + [ MB t + [ 4>tdQ t + [ vA-Stdt, 
Jo Jo Jo Jo 

where the last term represents the amount of money received as dividends. The cost 

process C t associated to the strategy is just the difference between the value of the 

portfolio IT and the revenues V t , C t = Ut — T> t , and the trading strategy is said to be 

self-financing if C t = 0, because in that case there is no net cash flow entering or leaving 

the replicating portfolio [30]. When C t = the variation in the value of the portfolio 

dU t coincides with the change in the capital gain dV t : 

dIT = u t dS t + <j) t dBt + ifitdQt + vA-Stdt. (9) 
Now we have to demand to that SDE (Q matches ([7]): 

dtPt-dt + dsPt-dS t + ^S 2 t d 2 ss P t -dt + APt-dh = 
u t dS t + MB t + iptdQt + vA-Stdt. (10) 
In an analogous way to dPt, dQt follows a SDE, 

dQ t = dtQt-dt + dsQt-dSt + ^alS 2 d 2 ss Q t -dt + AQ t -dI t , (11) 

where 



AQ t = Q(t, S t , 1) - Q(t, S t , 0) = Q b (t, S t ) - Q a (t, S t ). 
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d t P t - 



1 



a z t _S 2 t d 2 ss P t . - u t 8 t St - r<j> t B t 



dt 



(u t - d s P t - + ilhdsQt-) dS t - (AP_ - V*AQ t _) dl t . (12) 

In order to transform (Tl2l from a SDE into a deterministic partial differential equation 
we must guarantee that terms containing the stochastic magnitudes dS t and dl t cancel 
out, by selecting the self-financing trading strategy. Therefore, we must demand that 

u t = d s P t - - iptdsQt-, 
a condition named delta hedging, and also that 
AP t _ 



which is usually referred as vega hedging. The previous hedging conditions reduce (112 
to 

dtPt- + \a 2 t _S 2 d 2 s P t ^ - 5 t ^S t d s P t - = 



rd t B t + 



2 

AP t _ 



dtQt- + -a 2 _S 2 d 2 ss Q t - - St-StdsQt- 



(13) 



AQ t _ 

an expression that still involves (/>tB t . We can make it disappear just by using the 
definition of the portfolio in (jSJ), the predictable character of 4> t , the delta hedging and 
the vega hedging, altogether, 



j t B t = P t _ - v t S t - 4>tQt- 



P_ 



d s Pt- 



Afj 
AQt- 



-d s Qt- S-, 



aq 4 _ 



The replacement of 0t-Bt in ffl3|) leads to 



d t P t . + -a 2 _S 2 «9§ 5 P_ - rP t _ + (r - 5 t .)3 t d s P t - 



AP 

Ag t 



at_Sfdi, s Q t - - rQ t _ + (r - 6 t -)S t d s Qt 



This formula implies the existence of a generic predictable process S t 
[T4"] . which uncouples the problem of finding P t and Q t : 

«9 t P_ + l -a 2 _S 2 d 2 ss P t ^ - rP t _ + (r - S t .)3 t d s P t . + H t AP t _ = 0, (14) 

and proves that the option Q t removes the risk associated with I t . Formula (fT4l is valid 
for the secondary option as well, merely by replacing P t with Q t . Note that, whereas (j3j) 
and affect this expression, it contains no reference either on jl t or A. Before discussing 
the financial interpretation of S t , let us briefly analyse formula (I14p . As before we will 
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use the notation E a (t,S t ) = E(t,S t ,0) and E b (t,S t ) = E(t,S t ,l). When t < r the 
differential equation (TH1) reduces to 



d t P a + \<J 2 a S 2 d 2 ss P a - rP a + (r - 5 a )S t d s P a + H Q (P b - P a ) = 0, (f 5) 

whereas for t > r we will have 

<9 t P& + \°lS 2 t d 2 ss P b - rP b + (r - 8 b )S t d s P b + E b (P b - P a ) = 0. 

Note however that in this case P b must fulfil the classical Black-Scholes-Merton equation 
with continuous-time dividends, i.e.: 

8 t P b + ^a 2 b S 2 d 2 ss P b - rP b + (r - 8 b )S t d s P b = 0, (16) 

therefore we have a first constraint for E t , E b (t,St) = 0. Finally, since Pt is a right- 
continuous stochastic process, we have that P(r, 5 T , I T ) = lim^ r Pb(t, St). 

We will concentrate now on S 4 . To this end, let us introduce another portfolio 



^ t = u t S t + t 5 t + P t , 

based upon the following self-financing trading strategy 

u t = - d s P t -, 

~ S t d s Pt- - Pt- 

0i = % ' 

On the one hand we have that — 0, and on the other hand 

d^ t = D t dS t + MB t + dP t + ZtSt-Stdt 

= - d s P t -dS t + [(r - 8 t _)S t d s P t - - rP t _] dt + dP t 

dt + AP+dL 



d t P t „ + l -a 2 t _S 2 t d 2 ss P t - - rP t _ + (r - 8 t „)S t d s P t - 



= AP t _ {dl t - E t dt) , 

where we have used ([7]) and (fi4|) . Let us analyse now the right-hand side of the last 
equality in the previous equation for t ^ r, and assume that for some range of values 
of its arguments E a (t,S t ) < 0, and therefore E t < 0. Moreover, let us restrict our 
considerations to a subset in which we have simultaneously E t < and AP t _ > — note 
that both quantities are predictable. The only thing we ought to do in order to obtain 
what is known as a "free lunch" , i. e. sure earnings without capital exposure, is to wait 
until convenient market conditions are reached and costless construct ty t - immediately, 
because then d^ t > 0. The possibility of performing this kind of investment is named 
arbitrage opportunity , and it constitutes an undesirable feature of any market model. 
Even if E a = we would get a "free ticket" for a lottery, d^t > with probability 
Xdt — the regime change either takes place — and d\I/ t = with probability 1 — Xdt 
— or not. Note that if AP t _ < the trick is to compose the reverse portfolio. Finally, 
the conclusion is valid for AP t _ = as well, because P(t,S t ,It) and E(t, S t , It-) are 
continuous in t and S t . E a (t,S t ) > is the only strong constraint on this magnitude, 
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which depends in essence on how every investor measures volatility and dividend risks, 
that are not traded in the market. 

Since S a cannot be unambiguously settled, the market is not complete strictly 
speaking. Therefore one must resort to some sound financial criteria to choose H a . 
A conventional choice is the one that removes the so-called statistical arbitrage that 
appears when the growth of the expected value of a null portfolio is different from zero. 
In our case this is attained under the measure P if S a = A, since then for f^rwe have 

E p [d*t| T t -\ = (A - E a ) AP 4 _dt = 0. 
4. Perpetual American vanilla options 

The main assumption that determines the range of validity of equation (fl4j) is the 
irrelevance of past values of the underlying in the valuation of the derivative. This 
constraint is satisfied for any option whose value at the exercise time, t*, depends only 
on the present price of the underlying, S t * , and t* itself: 

P(t*,S t »,I t »)=X(S t *)h^ T . 

When the option can be exercised at the end of the contract lifetime only, t* = T, the 
exercise time is deterministic, and the option is said to be European. If the option 
can be exercised at any time before expiration it is called American, and t* becomes a 
stochastic magnitude as well. Note that the contract is always worthless after maturity, 
t* > T, and therefore the option buyer must decide under which conditions the option 
can be optimally exercised before this deadline. Since the decision may be based on S t 
and the natural time variable, the time to maturity, T — t, one deduces that the optimal 
exercise boundary, the stock price for which it is better to exercise the option than to 
keep it alive, will be of the form H t = H(t,I t ), with the usual auxiliary expressions 
H a {t) = H(t, 0) and Hb(t) = H(t, 1) assumed. The very definition of H t entails that for 
any live American option it holds that 

P(t,S t ,I t )>X(S t ), (17) 

and that H t must be settled in such a way that relationship 

P(t,H t ,I t ) = X(H t ). 

is satisfied. As a consequence, the computation of the option price and the optimal 
exercise boundary must be tackled contemporarily. Here lies the root of the difficulty in 
pricing American derivatives: first one must solve the partial differential equation (fl4"I) 
for any possible boundary condition, P(t, S t , It ; G t ), 

P(t,G t ,I t ;G t )=X(G t ), 

and then infer the value of H t from some other financial argument. Note that option 
holders try to delay the exercise of the option as much as they can, since the worth of a 
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live option is never lower than the pay-off of the contingent claim, (1171) . Then H t is the 
function that maximizes the price of the option with respect to the boundary condition: 

P(t,S t ,I t ;H t )^P(t,St,I t ;G t ). (18) 

When the pay-off function is regular enough, the above demand is equivalent to the 
classical smooth pasting condition that states J2j H] : 

dX{S) 



d s P(t,S,I t ;H t 



(19) 

S=H t 



is=m ~ dS 

In conclusion, the option price must be continuous with continuous derivative in the 
asset price when it crosses the boundary. 

Few closed solutions for this general problem are known even when the stock price 
follows a standard geometric Brownian motion. In particular the question stands open 
for the most ubiquitous pay-off function, 

X ± {S) = max(±(S-K),0). (20) 

This kind of derivatives are known as vanilla options, and the plus (respectively 
minus) superscript stands for call (respectively put) options. The difficulty deceases 
dramatically if one can deduce the way in which the optimal exercise boundary depends 
on time. This is the case when the option is perpetual, when it does not mature. 
It can be objected that perpetual American options have limited practical interest 
since actual derivatives commonly expire. However, since they can be understood as 
an approximation of a far-from- maturity contract, perpetual options may help in the 
valuation process when theoretical prices cannot be computed [10] . Moreover, as we will 
stress later, the emerging problem is still akin to physics because it involves mechanical 
or thermodynamic concepts like equilibrium and stability. 

We will translate now the absence-of-expiration assumption into our expressions 
by letting T — > oo. This will remove any explicit temporal dependence from P t , 
P(t,St,It) — > P {St,It) and from H t , H(t,I t ) — ► H{I t ). So, the remaining implicit 
time dependence in H t is fruit of a possible future change in the dynamics exclusively, 
and therefore 

Ht = H a l t<T + Hfjlt^r, 
with H a and Hj, non-negative constants to be determined. 

5. Main analytic results 

In this Section we will present the solution to equations f|T5|) and ffT6l) . with E a = A, 
when the pay-off function is ( 1201 . The detailed derivation of the expressions to appear 
can be found in|Appendix A|for the case of call options, and one may use the procedure 



developed there as a guideline for finding put formulas. 
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with 



ft 



ft-^ 



± m 



K. 



K) 



Here and hereafter we will denote by ft b the following combination of the model 
parameters r, 5 a ,b, and a a ^'- 

1 



ft, 



a 



a.b 



A,6 ± y/Pafi + 2r < b ) 



(21) 



where 9 ab is a shorthand for 

a ,b = r - 5 a , b - a 2 a b /2. 

Note that ft b ^ 0. In fact, one can prove that ft b ^ 1 and > (3~ b ^ —2r/a^ b , where 
equalities hold if and only if 5 a ^ = 0. 

The relative values of ft and j3 b play a key role in the evaluation of the price of 



the perpetual option before the regime change. When \(3^\ ^ \ft\, i-e. when (3^ -p (3 f 



< 



one finds that necessarily Hf ^ H b , and the value of P<t{S; Hf ) can be encoded in a 
single expression: 

P±(S ] H±) = ±(H±-K) 



±{Ht ~ K) 



A 



_S_ 



S 

Q± 



Hi 



S 



(s 



^ H a)i 



where we have introduced 7^ ^ which incorporates A: 
1 



it 



-0* ± y/0* + 2(r + \)al 



(22) 



(23) 



and constant which depends on the remaining parameters of the market model: 

r 1 



(ft - ft) 



ft 



(24) 



and cancels if and only if ft = ft- The value of the optimal exercise price cannot be 
obtained in a closed form valid for any value of the parameters. As we have explained 
above one can recover by demanding the so-called smooth pasting condition ffT9l) 
on the solution ff22l . This one is the approach we follow in the Appendix because it 
naturally provides ordinary differential equations with Cauchy boundary conditions on 
S = H^. Here we will show how (jl~9l) is just a consequence of the general maximal 
principle (|T8l) which in this case takes the form: 



0. 



(25) 



We will drop the subscript in St hereafter for notational brevity. 
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If one differentiates (1221) with respect to H^ gets 



g 



Hi 
Hi 



Hi 



where we have defined the auxiliary function Q (77), 

Q±(r)) = A^rf-ft - rf^ =F B ± , 
and the constant positive-definite coefficients [jj 

1 /? 



7^ 



5 



± 



7± 



^-1 
A 



7? 



(26) 



(27) 



Therefore we have to look for 77^, ^ ± (?7^ = ) = 0, because then the value of H^ will follow, 
Ht = VoHt 

It is shown in the Appendix that function G + (i]) = has one and only one zero in 
the interval < 77 ^ 1, but the same kind of analysis performed there leads to the 
conclusion that Q~(rf) = has a single solution in 1 ^ 77 < 00. Thus the problem is 
properly settled. 

When |/3^| < |/^| the optimal exercise price for vanilla calls (respectively puts) 
is higher (respectively lower) before the regime change than after it. This makes 
feasible that a change in the dynamics compels the holder to early exercise the option 
immediately because the sudden change in the optimal boundary value. This fact leads 
to the existence of two different pricing formulas depending on the present value of the 
underlying: 



± 



P^S;Ht] 

±_* *_ 

\ + r 7 ± _ 7 t 

la la la 



s a 



X + 5 a 



-m 



it 



7 



A + r 
t-1 



-K 



^6 - 1 A 



± 



It 



_s_ 

Ht 



1-7?)^ -7?) 



±\ (7a -7?) 



5 



it 



it 



(S § flf ), 



(2? 



and 



AiT 



± 



S a 



X + 5 a 



-Hf- 



X + r 



-K 



-Y 



7" 



"A + r( 7a ± -7?)(l-7j)(A 



S 

Ht 



H, 



±\ 7. 



5 

1? 



7a= 



A + 5 a 



A + r 



(29) 



|| One can consider different definitions for the auxiliary function and the associated coefficients. We 
kept the same notation appearing in 19J. 
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Both expressions match when S = Hf, in a very smooth way: one has to compute the 
fourth derivative of the option price with respect to the spot price before a discrepancy 
is found. The value of Hf remains elusive: the optimal exercise price solves an equation 
that in the general situation is transcendental again. In this case the use of the 
maximization criterion ( 1251) instead of the matching condition (|T9l) is very illustrative 
as well. Since ( 125]) is a global condition, it affects the whole solution. Then one may 
perform the computation either on fl28|) or on (1291) with the same result, namely, 



with 



ft \ + r\Ht 



S 



la 



Ht 
Ht 



Ht 



nHx) = x 



" 7? =f C 



-7? ± D ~- 



and 



C 



D z 



± 



it 



Pt 5 a 



1-jfft-lr 



± 



A 



r(ft-li)il-li) 



(30) 
(31) 



positive magnitudes. As before, there is a single solution to the equation T~C ± (xo) = 



that belongs to the interval 1 < x 
compute though 

H a = Xo H b ■ 



±i 



< oo. After obtaining the value of Xo we can 



Note finally that in expression (1281) it was assumed that ft ^ 7, 



When a numerical 



concordance in the parameters leads to ft = ft, we have merely to take the limit 
ft — > ft in the corresponding expression and the right solution is obtained: 



A 



± 



S a 



K 



it 



Ht 



'X + rft 

(it) 2 



7? 



it 



A + r 

■27a* 



-if 



g 

7± 



7a± 



+ 7^1n 



7± 



Ht_ 
Ht 



{la 



1 

-7?) 



it 

s 

Ht 



it 



S 



g 



7a± 



(32) 



V la a -It I V a / \ 

The question is that the case in which ft = ft does not constitute a critical point for 
the process, it only affects the functional form of the solution. There are some other 
situations with deeper financial relevance that deserve a closer study. We will introduce 
and illustrate them in the next Section. 



6. Discussion 



Let us analyse in the first place the consequences on call option prices of a stoppage 
in the payment of dividends, 5 a 7^ and 5b = 0. After the change we have ft = 1, 
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El 



oo, and P, + — > S 1 . This ensures that the right choice here is the solution for 



which H£ ^ H£ . So we have to write equation (1221) for the case in which Pi — 1: 



5 n 



-h: 



K 



where we have taken into account that here £ + 



S 

Hi 

6 n 



it 



+ 



A 



H]. Since /3 6 + 



S. 



(S < Hi), 
1 implies A + 



OO 



and P + —>■ oo, it is somewhat simpler to recompute (1251) directly, 



'la 



s 

Hi 



la 



'a 



1 5, 



it 



9H+ 

what leads to the following algebraic expression for 

ii(i + x/si 



A + 5 a 



-K. 



Note that this expression for Hi leads to Pi < S. This result has its origins in the 
uneven rights that option holders and stock owners have before the change due to the 
dividend payment. Remember that after the change P fc + = S, so the call holder must 
face the following dilemma: if the option is not exercised there is still a chance for 
the sudden revaluation of the option price, since Pi > Pi, but at the same time, the 
longer he/she waits the higher is the amount of dividends he/she will never receive. 
The optimal exercise price Hi comes from a trade-off between these two competing 
reasons. Thus the difference disappears if we consider the limit 5 a — > 0: we can see how 
we recover Hi — > oo, the option is never exercised also in this case, and Pi — > S as 
well. Therefore if the share pays no dividend at all, American call options quote as its 
underlying even if the volatility changes. 

Let us assume next that 5 a = but with 5b > 0, the stock can start suddenly to 
pay dividends. Since Pi > Pi the appropriate pricing expressions follows from (1251) 
and (1291) . and Hi / Hi should be a zero of the auxiliary function H + (x)- The first point 
that deserves attention is the fact that in this case we have C + = 0, D + > and, as 
a consequence, H + (x) is strictly positive. Then the only way in which one can satisfy 
the maximal condition (j23|) is by demanding that Hi — ► oo. This conclusion must hold 
regardless the value of the rest of the parameters in general, and Pi > 1 in particular. 
The validity of this result was already conjectured in [19J, but no prove was given there. 
Note that in this case the price of the call differs from the underlying, < S, due to 
the existence of the post-change scenario, P t + < S: 

a Pi 



and 



pi(s) 
+(Hi- 

Pa + (S) 



A + r 7 + 
A 



K) 



la (li 
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X + £ + \H 
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0t 
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-K 
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Hi 
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+S - K, {Hi < S < oo) 

X + r b 
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Only in the limit A — > 0, the solution P+ = S is recovered. Figure [T] shows how H+ 
increases when one considers lower and lower values for 5 a , until the limit 5 a = is 
attained. In this case the price never intersects the pay-off function. 




S/K 

Figure 1. Option prices for perpetual vanilla calls under dividend risk. We represent 
the price of the option, in terms of the moneyness (S/K), for different values of 5 a . 
We observe how both the optimal exercise price and the option price itself increase as 
S a tends to zero. When S a = the option price never touches the pay-off function, 
depicted in a solid black line. We have used the following values for the parameters: 
r = 4%, a a = cr fc = 10%, 5 b = 2.5% and A^ 1 = 2 years. 

Another aspect of the previous solutions that it is worthwhile devoting attention 
to is the behaviour of those expressions out of their range of validity. As we have 
pointed out in the previous Section, and partially shown in the Appendix, equation 
H ± ( y Xo) — has a unique solution in the interval 1 < x ±l < 00 when \(3^\ < The 
fact is that if one considers arbitrary relative values for f3^ b and expands the domain 
of H ± (x) to the whole positive real axis, it always presents a single zero, because 
\im x ^ TC ± (x) < 0, lim x ^oo TC ± (x) > 0, and the function is either monotonous or has 
a single local extreme. Therefore one can compute Hf = XqH^, H ± (xo) — 0, when 
\Pt\ ^ \Pb\i an d compare it with the right value = r/^H^. The same kind of 
analysis can be performed on the two associated pricing expressions: P^ 1 , the correct 
one in the present circumstances, found in (f22l) . and the heuristic one obtained by 
using formula (1251) . Note that solution (1291) does not apply in any case. This is not a 
mere academic exercise, because both expressions solve the same equation, namely (1A.6I) 
for call options, and still obeys to a maximization principle. The divergence emerges 
as a consequence of the different boundary conditions to be satisfied: in particular it is 
no longer true that P^(S = Hf) = ±(Hq — K). Since this condition is deduced from 
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arguments on optimality, one can consider the heuristic expression as a non-optimal 
solution. In Figure [2] we can observe the relative difference between pricing expressions, 
(P~ — P~)/P~, for a perpetual put option in a market model with /3~ = —8.0 and 
fit = —1.28. The heuristic solution underprices the put in a quantity that reduces for 




0.5 0.6 0.7 0.8 0.9 1.0 

S/K 



Figure 2. Relative price change for a perpetual vanilla put under volatility risk. We 
represent the price of the option, in terms of the moneyness (S/K) , for different values 
of A -1 , expressed in years. We analyse here the difference between expressions (|22|) 
and (|28p when a sudden increment in the volatility of the stock is presumed. We have 
used the following values for the parameters: r = 4%, 6 a = Sb = 0, a a = 10% and 
cr b = 25%. Note how the right solution dictates a higher price. 

large and small values of S and A, because then the mismatch in the boundary reduces 
its relevance. In any case, the error introduced by the heuristic pricing expression in 
the terms analysed is relatively small and bounded. In fact, as we show below with a 
different example, if we simply superimposed the two prices in a single plot it would 
seem that they just collapse. Besides P~ ^ P~ we have also H~ ^ H~. This is tenable 
because the sudden change in the heuristic solution at S — H~ , clearly noticeable in 
Figure [2j P~ also crosses the pay-off function at S > H~, what reinforces the heuristic 
character of the solution. 

Let us consider now the opposite situation in which one tries to use ( 1221) instead 
of expressions (!28!) - (!29l) when < |/3^ |. The immediate conclusion is that in this 
case ( 1221) does not solve the appropriate differential equation, because P b — which 
plays the role of an external interaction — is not properly settled, but right boundary 
conditions are to be satisfied eventually This may result into incompatibility. Note 
that lim^o Q ± (i]) ^ 0, and lim^oo Q ± {rf) ^ as well. This implies here that G ± (rj) has 
either two zeros or none. If ^ we have stated that the problem has always a 
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univocal solution, but when < \(3j^\ it is possible that no extremal solution can be 
defined. This problem was already detected and depicted in [19]. In Figure [3] we observe 
how, for a given market conditions, no heuristic solution can be found beyond a critical 
value of A, A m 0.5094. Once again the value of the heuristic optimal exercise price is 
smaller than H~ , and since the heuristic solution fulfils all the boundary conditions, it 
overprices the value of the option. 




S/K 

Figure 3. Option prices for a perpetual vanilla put under volatility risk. We represent 
the price of the option, in terms of the moneyness (S/K), for different values of A, 
expressed in years -1 , when one awaits a severe reduction in the volatility of the stock. 
We observe how right prices (solid lines) coincide optically with heuristic ones (dots) 
up to the threshold A « 0.5094 is attained, even for large stock prices — see the inset. 
For higher values of A we have only the correct solution. We have used the following 
values for the rest of parameters: r — 4%, S a = 5b = 1-75%, a a — 40% and Ob = 25%. 

Apart from the surprising accuracy of the heuristic solutions, the main conclusion of 
the previous discussion is that the right H^ 1 is extremal among different configurations, 
Figure HI Therefore one can optionally interpret that a solution cease to be valid because 
it conveys to an exercise price which is not as optimal as it could be, like stable and 
metastable equilibrium states in thermodynamics. 

Note finally that when (3^ = j3^, not only = but also P^ = P^, no matter 
the values that the market parameters take individually. 

7. Conclusions 

In this article we have obtained the complete set of analytic pricing expressions for 
perpetual American vanilla options within a market model that may present a single 
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Figure 4. Optimal exercise price difference. We represent H~ — H~ as a function 
of the "order parameter" |/3~ \/\f3^\, for different values of A, expressed in years -1 . 
We observe how the true value is always maximal, the difference is bounded, and H~ 
ends in a point where its derivative is infinite. In the confection of the plot we have 
varied a a and left unchanged the rest of parameters: r = 4%, S a = 8b = 1.75%, and 
<7 6 = 25%. 

change of regime in both the dividend-payment rate and the volatility level of the 
underlying stock. We further assume that the value of these parameters after the change 
are known quantities. The problem is relevant from the point of view of mathematical 
finance and shows interesting connections with physics. 

The most exciting financial result is perhaps that the price of a perpetual American 
vanilla call when the underlying does not pay dividends may differ from the spot price of 
the stock, regardless the fact that the option will be never early exercised as long as the 
market conditions remain fixed. The risk associated to the fact that this non-dividend 
policy can be eventually revised in the future reduces the value of the option. 

On the realms of physics one finds certain parallelisms between some properties 
shown by this model and phase transitions in thermodynamics: when one analyses the 
issue before the regime change, two different formulas are found for both the price of the 
option and the optimal exercise boundary, and the appropriate expression to be used is 
elucidated on the basis of the value of a some sort of "order parameter" . 

This analogy with thermodynamics extends to the mutual relationship between 
the two available exercise prices. Like state functions in thermodynamics, exercise 
prices characterize equilibrium states in the sense that they are solutions to extreme 
value problems. When two possible values for a function of state coexist, the one that 
is smaller (or larger, depending on the state quantity analysed) is seen as the stable 
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solution, whereas the alternative solution is named as metastable. In our case this kind 
of behaviour is observed, and one can trace the reason back to the same origin: each 
solution reflects different environmental conditions, and only one of them incorporates 
the actual constraints. And as in physical problems, the metastable solution may become 
eventually unstable and just disappear. We have shown how this phenomenon was on 
the basis of the spurious bizarre results reported in [19] . 

Finally, we have illustrated the major results enumerated above with several 
graphical examples corresponding to practical market situations. 
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Appendix A. Detailed computations 

Let us consider the case of call options in detail. The equation to be solved if the regime 
has changed is ([TBI) : 

^ 2 ^ + (r-W^-rP+ = 0, 

whose general solution reads 

P b + (S;H+) = K 1 S> 3 > +k SK, 

where f3^ first appeared in (12T|) . The applicable boundary conditions are three in this 
case: 

limP+ =0, (A.l) 
lim P+ =H+-K, (A.2) 

c\P + 

lim -A. = 1. (A.3) 



Since ^ ^ we have that k = by virtue of (lA.ip . Ki is obtained after flA.2[) . 
« a = (H+ - K) 



and therefore HU1: 



S 



K = {Ht-K)( — \ , (S^H+). (A.4) 



The value of H b emerges from (1A.3I) : 



Ht = ^K. (A.5) 
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It can be proven that {3^ ^ 1, and that (3^ = 1 is equivalent to 5^ = 0. Therefore we 
will have K < < oo for p£ > 1, and for (3^ = 1 the value of if^ diverges: there is 
no optimal boundary, the option is never exercised and quotes as the underlying stock, 

lim P+ = S. 

Note that (1A.4j) is valid for S ^ H^~, whereas the value of the call for 5* > is just 
S — K. Some remarks about this statement can become very revealing for a better 
understanding of future developments. The exercise price of a call marks the threshold 
above which it is better to execute the option rather than to keep it alive. However, this 
does not ensure that one will be able to exercise the option exactly when S = H^~. Let 
us think, for instance, in the case of a market with quantified quotes, and observe how 
in general the executed price will exceed the optimal exercise price unless it coincides 
with one of the admissible prices. In our case, the process S t is continuous but the 
optimal exercise price is not, it changes from if+ to H^~, what can end in the same 
result if H+ > H+. 

Let us consider now the call price before the change, the one we are primarily 
interested in. The differential equation is 
1 c\ 2 P + c\P + 

2^ 2 ~o# + (r " 6a)s ^r ~ {r+ x)Pa+ = _An+ ' {s ^ H « l 

and therefore one must choose the right expression for P b + . Two well different scenarios 
emerges, depending on the relative values of and . If ^ H£ we will have a 
single equation to solve 

r 2aS2 ^ +{r ~ 6a)s ^~ {r+x)Pa+ = - x w- R ) >( A - 6 ) 

for all S ^ H£, whereas when H£ > H^~, a set of two different equations appears, 
equation (1A.6I) itself for S ^ H^~, and 

1 d 2 P+ dP+ 

2 ff ^ 2 ^# + (r " 6a)S ^S~ ~ {r + X)Pa+ = ~ HS " K) > (A ' 7) 
for H£ < S ^ H*. The point is that the value of is obtained after the resolution 
of the problem. Therefore we will start our analysis by assuming that H£ ^ H£ . The 
general solution of (1A.6|) is 



P+ = k 3 S^ +Ka ST' +YT¥ (Ht-K) ^ j ,(S< K ^ H+), (A.8) 

where 7^ are defined in (1231) . £ + in fl2M . and we are assuming for the moment that 
7a" 7^ Pb '• Now we can impose the boundary conditions 

lim P+ = 0, (A.9) 

lim P+ =H+-K, (A. 10) 

S-+H+ 



In fact, most of real markets behave in that way, because they demand a minimum variation to 
change a price. This tick size may be as small as a cent, but this does not invalidate the argument. 
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(A.ll) 



From ( 1A.9I) we deduce k 2 = 0, and ( lA.lOj) leads to 



«3 



and then 1T91: 



A 



A 



A + £+ 



HI 



la 



S 



la 



s 



Ht) 



S 



la 



,(S^tf+).(A.12) 



The value of H£ must fulfil the transcendental equation that follows from the smoothness 
condition flA.111) . 



7+ - 1 " 



2A 



<ft(ft-% 



K 
H+ 



0. 



where H£ was introduced in (IA.5I) . When / 1 we can define the auxiliary function 



g+( v ) = A + r] 



Br 



7] 



with A + and B + bounded and positive constants — see (1251) and f[2"T|) in the main text 
for their definitions. The analysis of this function is relevant in our framework because 
fulfils Q + (H^ / H^r) = 0. Since we have assumed that ^ H^~, we are interested 
in studying the properties of Q + {rf) in the range < 77 ^ 1. The values of the function 
in the interval endings are lim^o {? + (??) — — °o and Q + {1) = £ + B + /X. G + {f]) shows a 
single extremal value located outside the segment r\ G (0, 1]: 



Vm 



ll 



> 1, 



-ft (VM)-ti- 2 < 0. 



7+-1 

which is a maximum since 
r d 2 £+ 

The above properties compel equation Q + {j]q) = to have at the most one solution 
for which < T)q ^ 1. The necessary and sufficient condition for the existence of such 
solution is 

G + (l)^0^£ + >0^ft>ft. 

Therefore, condition ^ grants that /?+ ^ and, incidentally, that 7+ > is 
satisfied in (TA~T2l . 

In fact, we can show how > H£ leads to /?+ < /3^" as well. Let us assume 
that if+ > if^: in this case we have to solve (1A.6I) and (1A.7I) . The general solutions if 
7a ^ ft are: 



P+ = K 5 S^ + K 4 S 



la 



A 

A + £+ 



x + - K) 



S 



ft 



b )i 
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P+ = K 7 S^ + K 6 S la 
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-S 
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X + 5 a A + r 
with the five boundary conditions to be satisfied listed below: 



K, (H+<S^H+), 



After imposing constraints (lA.13l) - (1A.16j) one gets 
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lim P+ 


= 0, 


(A.13) 


lim P+ 


= lim P+, 


(A. 14) 


dP + 
lim + ^Tf 


dP + 

S|# 6 + db 


(A.15) 


lim P+ 


= K - K, 


(A.16) 


dP+ 
d6 


= 1. 


(A.17) 



A + 5 a A + r 
Condition (1A.17|) leads to a new transcendental equation: 



k - 1; 



(l- 7 -)(/3+- 7 -) 



a * +r— + -^^5, 



„ - 0, 



5 



h + 1 - 7- 

a 'a 

and one can define a second auxiliary function in order to analyse the problem: 

since H + (H+/H^) = 0. As before C + and D + are bounded and positive constants, 
defined in the main text, equations ( l30l) and ( J3TT) . The values at the extremes of the 
interval of interest are limine 7i + {\) = —00 and 



la 



*5« W = w- 7 - )W -i )r 



>o, 
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provided that /3+ < f3^ . When 5 a ^ function 7i + (x) nas a maximum at xm, 



Xm 



since 



- 1 r ' 



d 2 H+ 

x^xm dx 2 



lim = 7" ( X m) 7 " 2 < 0, 



whereas 7i + {x) is a decreasing function if S a = 0. In any case, one can conclude that 
7i + (x) has a single zero xt irrespective of the maximum location which, if 

exists, may be placed either inside or outside this region. 

The resolution scheme for perpetual American vanilla puts is very similar. We have 
decided not to include it here for brevity reasons. 
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